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1.  Introduction 

In  a  previous  report  [l] ,  we  presented  the  results  of  our  investiga¬ 
tion  up  to  June,  1964  of  waves  incident  upon  a  beach.  Particular  emphasis  was 
placed  on  the  wave  height  at  the  undisturbed  shore -line  and  the  motion  of  the 
shoreline,  including  the  maximum  run-up  distance.  For  time  harmonic 

waves  of  small  amplitude  incident  from  the  open  sea,  formulas  for  these  quantities 
were  derived  based  upon  the  linear  theory  of  water  waves.  The  nonlinear  theory 
was  used  to  find  the  range  of  validity  of  these  formulas.  In  addition  a  numerical 
method  was  devised  for  the  calculation  of  the  wave  motion  on  the  basis  of  the 
nonlinear  shallow  water  theory.  This  method  was  applied  to  waves  on  a  uniformly 
sloping  beach.  The  results  showed  fair  agreement  with  the  analytical  results 
for  waves  of  low  frequency  but  not  for  higher  frequencies. 

From  June  1964  to  June  1965  we  have  improved  the  numerical  method  by 
using  a  finite  difference  scheme  of  higher  order  accuracy  described  in  section 
3  of  this  report.  With  the  improved  method  we  have  made  extensive  calculations 
which  agree  much  better  with  the  analytical  results  and  which  do  so  up  to  higher 
frequencies.  However  the  discrepancies  between  the  numerical  and  analytical 
results  still  occur  at  frequencies  somewhat  below  those  at  which  the  analytical 
results  cease  to  be  valid.  Although  we  have  not  found  the  cause  of  these  dis¬ 
crepancies,  we  can  conclude  that  at  the  lower  frequencies  the  numerical  results 
are  reliable.  Some  of  the  results  are  described  in  section  5  of  this  report. 

We  have  also  developed  the  analytical  theory  of  wave  breaking  and  bore 
formation  and  growth  on  a  uniformly  sloping  beach.  This  theory  is  based  upon  the 
nonlinear  shallow  water  theory.  By  employing  the  method  of  characteristics  in 
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section  6  we  determine  how  an  incident  wave  becomes  distorted  and  finally  reaches 
the  breaking  point.  Then  in  section  7  we  utilize  the  theory  of  weak  shock  waves 
to  determine  how  the  bore  forms  and  grows.  This  analysis,  being  based  on  the 
hypothesis  that  the  wave  is  vaak,  cannot  yield  all  the  properties  of  the  bore, 
especially  those  at  the  shoreline.  Nevertheless  it  does  determine  how  the  height 
of  the  bore  increases  rapidly  at  first  and  then  more  gradually.  It  also  determines 
the  position  and  speed  of  the  bore  as  functions  of  time. 

This  theory  of  bore  formation  and  growth  on  a  sloping  beach  should 
form  the  basis  of  further  theoretical  and  numerical  investigations  of  water 
wave  run-up. 
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2.  Review  of  analytical  results 

In  order  to  compare  our  present  numerical  results  with  our  previous 
analytical  ones,  ve  shall  first  review  the  analytical  results  contained  in  our 
previous  report.  It  is  convenient  to  present  those  results  in  terms  of  the 
dimensionless  variables  which  we  have  used  in  the  numerical  work.  These  dimension¬ 
less  variables  are  the  time,  t,  the  horizontal  distance  x  measured  from  the  un¬ 
disturbed  shoreline,  the  undisturbed  water  depth  H(x),  the  disturbed  water  depth 
h(x,t),  the  shoreline  position  £(t)  and  the  horizontal  water  velocity  u(x,t). 

They  are  related  to  the  corresponding  dimensional  quantities  t,  x,  H,  h,  £  and  u 
by  the  scaling  relations 

x  =  x  D  cot  a,  h  -  h  D,  H  =  H  D 
t  =■  tfD/g)1/2  cot  a,  u  =  u(gD)1/2,  £  =  £  D  cot  a 

Here  a  is  the  beach  angle  at  the  undisturbed  shoreline  x  =  0  and  D  is  a  vertical 
length  scale  factor.  It  is  to  be  noted  that  the  horizontal  length  scale  factor 
is  D  cot  a.  See  Figures  la  and  lb. 

We  have  found  it  conven  lent,  to  choose  D  so  that  the  dimensionless 
undisturbed  water  depth  H  is  unity  at  x  =  -1.  In  addition  the  undisturbed 
depth  is  zero  at  the  undisturbed  shoreline.  Thus  H(x)  must  satisfy  the  two 
conditions 

2.  H(0)  =  0,  H(  — 1)  =  1 

It  is  assumed  that  the  water  lies  to  the  left  of  the  undisturbed  shoreline  in 
the  region  x  <  0  before  the  disturbance  arrives.  For  a  uniformly  sloping  bottom 
profile,  the  dimensionless  undisturbed  depth  is  simply 

3.  H(x)  =  -x. 
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Let  us  assume  that  the  incident  wave  height  [h(x,t)-H(x)Jinc  in  the 
deep  water  far  from  shore,  i.e.,  in  the  region  x  «  -1,  is  given  by  a  periodic 
sinusoidal  wave 

[h(x,t)-H(x)]lnc  =  a  cos(wt-px) 

Here  the  dimensionless  angular  frequency  w,  the  dimensionless  wave  number  p  and 
the  dimensionless  amplitude  a  occur.  They  are  related  to  the  dimensional  angular 
freqiency  u>  and  amplitude  a  by  the  relations 

5.  a  =  ^(D/g)1/2  cot  a,  a  =  a/D,  p  =  w2g -1D  cot  a  =  tan  a. 

The  amplification  factor  A  is  defined  as  the  ratio  of  the  maximum  water  depth  at 
the  undisturbed  shoreline  to  the  amplitude  of  the  incident  wave  far  from  shore. 
Thus 

6.  A  =  max|h(0,t)/a| 

t 

We  have  previously  determined  A  for  water  of  slowly  varying  depth.  This 
means  water  for  which  the  fractional  change  of  depth  in  a  wavelength  is  small 
except  near  the  shoreline.  When  this  condition  is  satisfied  and  lim  H(x)  =  1, 

X  — >  —  oo 

we  find  that  A  depends  only  upon  w  and  a  and  is  given  by 


2  2 

Here  z  ,  a  function  of  to  tan  a,  is  the  positive  loot  of  the  equation 

2  2 

8.  z  tanh  z  -  u>  tan  a 

In  this  case  it  is  not  required  that  H(-l)  =  1  since  D  is  chosen  as  the  undis¬ 
turbed  depth  at  x  =  -oo.  Graphs  of  A  as  a  function  of  u>  for  various  beach  angles 
a,  based  upon  (7),  are  given  in  Figures  3  and  7.  The  limit  of  A(w,a)  as  a  tends 
to  zero  is  easily  found  to  be 


< 
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9.  A(w,0)  =  v/^rtw 

Since  the  result  (9)  applies  when  w  is  fixed,  it  follows  from  (5)  that  ”  tends 
to  zero  as  a  does.  Therefore  (9)  is  the  asymptotic  form  of  (7)  for  low  frequencies 

w. 

The  result  (7)  was  obtained  in  two  ways.  First  it  was  deduced  from  the 
linear  theory  of  water  waves  in  water  of  any  depth.  Rien  it  was  rederived  by 
matching  the  result  of  the  linear  theory  in  deep  water  to  the  result  of  the  non¬ 
linear  shallow  water  theory  employed  near  the  shore.  The  latter  derivation 
also  yields  a  limitation  on  the  incident  amplitude  for  the  result  to  be  valid. 

It  implies  that  the  amplified  wave  height  be  small  compared  to  the  wavelength. 
Analytically  the  condition  is 


10. 


a  <  a 
-  max 


w  A(w,a) 


Here  A  is  given  by  (7).  Graphs  of  this  limiting  amplitude  as  a  function  of  u>  are 
shown  in  Figure  4  for  various  values  of  a. 

The  result  (9)  was  also  derived  in  another  way.  It  was  shown  to  be 
the  high  frequency  asymptotic  form  of  the  amplification  given  by  the  linear  shallow 
water  theory  in  water  of  arbitrarily  varying  depth.  Thus  the  result  (7),  which  is 
valid  for  high  frequency  waves,  has  a  low  frequency  asymptotic  form  which  is  the 
same  as  the  high  frequency  asymptotic  form  of  the  result  of  the  shallow  water 
theory,  a  theory  which  is  valid  for  low  frequencies.  Thie  shows  that  the  two 
results  are  asymptotically  equal  in  some  intermediate  frequency  rarge,  as  one 
should  expect  and  that  the  two  results  supplement  each  other. 

For  very  low  frequencies,  the  linear  shallow  water  theory  yields 
A  =  2  +  2w^b  + 
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Here  b  is  defined  by 


12. 


00 


b  = 


J  x[H‘1(x)-H‘1(»)]dx 
0 


In  the  special  case  of  a  piecewise  linear  bottom  profile,  the  equations 
of  the  linear  shallow  water  theory  can  be  solved  explicitly  for  all  frequencies. 
We  have  considered  the  special  case  in  which 
13.  H(x)  =  -x,  -1  <  x  <  0 

=1,  x  <  -1 

In  this  case  the  amplification  is  given  by 


14.  AKa)  +  J?  > 


1/2 


V 


Here  J  and  are  Bessel  functions  of  orders  zero  and  one  respectively.  This 
o  1 

function  agrees  with  ( II)  for  small  values  of  u>  and  with  (9)  for  large  values  of  w. 
A  graph  of  A  based  on  (l4)  is  shown  in  Figure  3  for  a  =  0.  It  is  clear  from  (14) 
that  when  a  is  small,  A  is  practically  independent  of  a. 

On  the  basis  of  all  of  the  preceding  results,  we  can  formulate  an  over¬ 
all  description  of  the  amplification  as  a  function  of  frequency  and  beach  angle. 

At  zero  frequency  the  amplification  is  two  and  then  it  increases  quadratically 
with  the  frequency  as  is  shown  in  (ll).  It  continues  to  increase  (see  Figure  3) 
in  a  undulatory  manner  as  is  shown  by  (14),  with  an  undulation  frequency  determined 

by  the  horizontal  scale  of  the  bottom  profile.  It  undulates  about  the  curve 

l/2 

A  =  (4jtw)  '  ,  which  it  approaches  as  u  increases.  Then  the  amplification  con¬ 


tinues  increasing  at  the  slower  rate  given  by  (7).  This  latter  transition 
shows  that  as  the  water  depth  becomes  appreciable  compared  to  the  wavelength, 


it  has  the  effect  of  diminishing  the  amplification. 

The  gradual  increase  of  amplification  with  frequency  just  described 
applies  only  until  bore  formation  occurs.  The  amplitude  at  which  this  occurs 
for  any  given  frequency  is  given  by  (10). 

A  cuantitative  theory  of  amplification  and  run-up  when  bores  form  has 
not  yet  been  developed. 

3.  Formulation  of  the  general  problem 

The  finite  difference  calculations  are  based  on  the  same  nonlinear 
shallow  water  theory  formulation  presented  in  [l] ,  which  we  repeat  here  for 
completeness. 

In  terms  of  the  variables  (l.l)  the  mass  conservation  equation  is 

1.  \  +  (uh)x  =  0, 

and  the  linear  momentum  conservation  equation  is 

2.  ut  +  uux  =  (H-h)x  . 

If  the  location  of  the  shoreline  at  time  t  is  denoted  by  x  =  £(t),  then  since 
the  water  depth  there  is  always  zero  we  have 

3.  h[£(t),tj  =  0. 

In  addition  the  velocity  of  the  shoreline  must  be  equal  to  the  water  velocity 
at  the  shoreline  so 

»•  ^  ■  »[«*>.»]. 

Equations  (3)  and  (4)  are  boundary  conditions  which  must  be  satisfied  on  the 


moving  shoreline 
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We  shall  consider  the  motion  only  in  the  interval  -1  <  x  <  £(t)  for 
t  >  0  since  we  assume  that  known  (subsonic)  waves  are  incident  from  x  <  -1 
for  all  t  >  0.  If  the  local  sound  speed  is  denoted  by  c  =  >/h  then  (l)  and 
(2)  can  be  written  in  characteristic  form  as 

5.  (u  +  2e)t  +  (u  +  c)(u  +  2c)x  =  Hx 

6.  (u  -  2c)t  +  (u  -  c)(u  -  2c)x  =  Hx. 

Equation  (5)  determines  (u+2:)  along  the  ’’positive”  characteristics,  dx/dt  =  u+c, 
which,  as  t  increases,  enter  x  >  -1  from  x  <  -1.  If  the  linear  shallow  water 
theory  were  valid  for  x  <  -1  and  H(x)  =  const  =  1  there,  then  the  general  solution 
would  be 

u  =  F(t-x)  +  f (t+x) , 
h  -1  =  F(t-x)  -  f(t+x). 

In  the  linear  theory  we  use  c  =  </h  =  1  +  (h-l)/2  and  so 
u(x,t)  +  2c(x,t)  =  2[l4F(t-x)]. 

Clearly  the  function  F  in  this  equation  represents  the  incident  part  of  the  linearized 
wave  motion.  Thus  we  take  for  the  boundary  condition  to  be  imposed  at  x  =  -1: 

7.  u(-l,t)  +  2c(-l,t)  =  2[l  +  F(t-l)  J 

The  water  is  to  be  at  rest  initially  so  that 
0.  h(x,0)  =  0,  u(x,0)  =0  in  -1  <  x  <  £(0)  =  0. 

Our  problem  is  thus:  to  solve  (l)  and  (2)  for  u(x,t)  and  h(x,t)  in  -1  <  x  <  £(t) 
for  t  >  0,  subject  to  the  conditions  (3)  and  (4)  on  the  unknown  boundary  x  =  £(t), 
the  condition  (7)  on  x  =  -I  and  conditions  (8)  at  t  =  C,  The  shoreline  position 
£(t)  must  be  determined  simultaneously.  We  shall  study  only  motions  in  which  bores 


are  not  present. 
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4 , _ Finite  difference  scheme. 

The  numerical  method  which  we  employ  is  of  higher  order  accuracy  then 
the  one  previously  employed  in  [l]  and  is  based  on  the  Lax-Wendroff  scheme  [2] . 
Again  we  use  a  uniform  spatial  net,  =  JAx  and  possibly  non-uniform  time  net, 
^k+±  =  +  *^ie  comPu"te(i  shoreline  position  at  time  t^  is  denoted  by 

=  and  ^en  xs(k)  is  define(i  as  the  net  V^int  for  which 

xs(k)  +  f  <  k  Z  xs(k)  +  * 

The  calculations  at  time  t.  are  then  divided  into  four  kinds  for  four  different 

k+1 

types  of  net  points,  as  follows: 

I)  interior  points:  <  x^  < 

II)  shoreline  interior  points:  xA  =  xg^ 

III)  incident  wave:  xq, 

IV )  shoreline :  £k+1 . 

At  points  cf  type  I  and  II  (see  Figures  8a ,b)  we  use  high  order 
accuracy  difference  approximations  of  the  equations  of  motion  ( 3. 1) ,  (3*2). 

Rather  than  write  down  the  messy  difference  equations,  which  are  not  instructive, 
we  shall  indicate  the  expansions  from  which  they  are  easily  deduced:  Thus  at 
any  point  P  of  Figure  8a  or  8b  an  application  of  Taylor's  theorem  yields 


h(P)  =h(P')  +Atht(P')  +^-htt(P')  +0(At3), 


u(P)  =u(P')  +Atut(P' )  +  “-uu(P')  +  0(At;5). 
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Using  (3.1)  and  (3.2)  ve  can  eliminate  and  on  the  right  hand  sides  above. 
Also,  differentiating  (3.1)  and  (3*2)  with  respect  to  t,  ve  can  eliminate  u. 

X  v 

and  h^t  to  get  finally 


2.  h(P)  =  h(P* )  -  At  ^  (uh) 


1  + 

a  [  aj 

Ip-  2 

¥[_  c 

*  +h  &  J 

0(&t3) 


3.  u(P)  =  u(P> )  -  At  |j(|  u24h-H)  + 


+  u 


S( |  u24h-H) 


"ST 


r 

-  p« 


+  O(At^) 


We  now  use  u,  h  and  H  at  the  points  L* ,  P*  and  R*  to  approximate  the  x-derivatlves 
appearing  on  the  right  hand  sides  of  (2)  and  (3).  Since  the  spacing  is  uniform 

p 

for  points  of  type  I  these  approximations  are  easily  obtained  to  accuracy  0(Ax“) 
by  using  centered  difference  quotients.  For  points  of  type  II  the  spacing  is  in 
general  non-uniform.  However  the  terms  multiplying  At  are  first  derivatives 
and  so  can  again  be  approximated  to  0(Ax  )  using  data  at  three  points.  The 
remaining  terms  are  approximated  to  within  O(Ax)  but  since  they  are  multiplied 

p 

by  At  ve  retain  third  order  accuracy.  One  should  have  no  difficulty  in  deriving 
the  difference  equations  from  the  above  description. 

The  difference  equations  at  points  of  type  III,  see  figure  8c,  are 
identical  with  those  used  previously  (see  eqs.  (37)  and  (38)  of  [l]).  They  are 
difference  approximations  of  the  characteristic  equation  (3*6)  centered  at  the 
midpoint  of  the  mesh  rectangle  in  figure  8c.  Combined  with  (3-7)  these  equations 
yield  u(P)  and  h(P)  to  third  order  accuracy. 

To  fit  in  the  position  of  the  shoreline  to  higher  order  accuracy  we 


again  note  that 
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At 


=  i(t)  +  At  *(t)  +  i(t)  +  O(At^) 

However  from  differentiating  (3.4)  and  using  (3.2)  at  x  =  £(t)  we  have 


k.  £(t)  =  [H(x)  -  h(x,t)]. 

so  that,  recalling  (3«^)> 


x=4(t) 


5.  t(t+at)  =  £(t)  +  At  u($(t),t)  +  ^-[H  (e(t)  -  h  (e(t),t)]+  0(At3) 


The  x-derivatives  in  (5)  can  be  approximated,  using  H  and  h  at  points  P*  and  R' 
of  figure  8b,  to  accuracy  O(Ax)  and  so  if  u(£(t),t)  is  computed  to  accuracy 

p 

0(Ax  )  we  get  £(t-»At)  =  £(R)  to  third  order  accuracy.  To  compute  u(£(t+At)  ,t4At) 
we  use  (3.4)  and  (4)  to  get: 

6.  u(£(t-*At),t+At)  =u(|(t),t)  +  At[Hx(^(t))-hx(^(t),t)]  +  0(At2) 

Thus  with  the  same  approximations  to  x-derivatives  used  in  (5)  we  obtain  u(R), 
see  figure  8b,  to  second  order  accuracy. 

The  stability  condition  and  interpolation  procedures  (for  >  xx+1  +  —) 
of  [1]  are  also  employed  but  we  do  not  repeat  them  here. 


5.  Results  of  the  calculations 

In  all  of  the  calculations  we  used  the  uniformly  sloping  bottom  profile 
(2.3)  and  uniform  x-spacing,  Ax  =  l/200.  The  incident  wave  function  used  in 
(3.7)  was 

1.  F(t-l)  =  a  sin  wt. 

Calculations  were  run  for  various  values  of  the  amplitude  r  and  frequency  u>  in 
the  following  ranges:  10  <  a  <  8  xlO  1  <  w  <  10.  Each  case  was  run  until 
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at  least  10  and  perhaps  as  many  as  14  full  period  vaves  had  reached  the  shoreline. 
In  most  cases  all  significant  transient  effects  were  absent  after  5  or  6  periods; 
however  in  some  cases  a  steady  regime  was  never  attained. 

ltypical  cases  in  which  time  harmonic  steady  states  were  established 
are  shown  in  figure  2.  In  figure  2  graphs  of  £(t)/a  are  shown  for  two  cases 
in  which  a  =  10  and  w  =  2  or  7.  The  steady  state  periods  in  these  cases  are 
Just  n  or  2rt/7.  After  about  the  first  1^  cycles  of  shoreline  motion  the  zeros 
of  £(t)  have  these  periods.  Significant  motion  of  the  shoreline  is  not  observed, 
in  either  case,  until  about  t  =  2.  However  this  is  just  the  time  required  for 
a  sonic  disturbance  to  propagate  from  x  =  -1  to  x  =  0  over  the  bottom  profile 
(2.3);  i-e. 


/  *  ‘I  M 


-1  -1 


-1/2 


dx  =  2. 


The  fact  that  the  signals  actually  arrive  slightly  sooner  may  be  attributed  to  the 
effects  of  the  nonlinearity,  near  the  shoreline,  which  increase  the  characteristic 
speed  (and  to  a  small  numerical  effect).  These  figures  also  show  that  the  maximum 
steady  state  runup  distance  may  be  exceeded  by  some  of  the  initial,  transient, 
runups. 

For  comparison  with  the  theoretical  results  we  define  the  computed 
amplification  as 


2. 


A  E 


T 


max 

<  t  <  T  +  — 

U) 


tit) 

a  9 


where  T  is  some  "large”  time.  If  this  quantity  varies  slightly  with  T  we  may 
use  an  average  over  several  periods.  In  Table  I  we  list  this  amplification  as 
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well  as  the  corresponding  minimum  shoreline  excursion  for  a  variety  of  cases. 

When  the  variation  of  these  quantities  over  the  last  three  computed  periods  was 
not  small  (say  less  than  5  ^> )  ve  list  the  range  of  observed  variation.  Thus, 
for  example,  the  case  u)  =  8  and  a  =  10’^  did  not  quite  reach  a  steady  state 
and  the  amplification  varied  between  12.9  <  A  <  14.2  in  the  last  few  periods. 

The  first  column  of  Table  I  is  plotted  in  figure  3  together  with 
theoretical  amplification  curves  for  a  =  0  and  a  =  tan"1  l/lO.  The  agreement 
is  very  good  up  to  w  =  6  (for  the  a  =  0  curve)  and  should  be  contrasted  with  the 
previous  calculations  (in  figure  5  of  [l]).  As  is  shown  in  this  figure  and  in 
the  table,  there  seems  to  be  some  kind  of  "resonance"  near  w  =  6.5  since  the 
amplification  becomes  very  large  there.  (There  are  lesser  oscillations  near 
w  =  4  and  u>  =  5-5)  Beyond  this  resonance  the  computed  amplification  values 
seem  to  lie  along  a  curve  considerably  above  the  theoretical  values  of  a  =  0. 

Of  course  as  the  frequency  increases  bores  will  finally  form  and  the 
present  calculations  must  fail  in  these  cases.  For  a  =  10  and  w  =  7  the  results 
show  quite  clearly  (see  figures  2b  and  5)  that  nothing  unusual  happens  and  in  fact 
that  bores  do  not  yet  form.  We  have  shown  in  figure  k  the  frequency  at  which 
the  calculations  seem  to  indicate  that  bores  form  for  various  values  of  incident 
wave  amplitude.  The  criterion  used  was  not  the  wild  oscillations  in  or  unsteady 
behavior  of  £(t)/a  but  rather  the  occurrence  of  negative  values  of  h  in  the 
calculations.  This  is  of  course  not  a  precise  determination  and  so  the  agreement 
with  the  theoretical  limit  of  equation  (2.10)  is  quite  gratifying. 

In  figure  5*  we  show  the  initial  wave,  for  u  =  2  and  a  =  10  ,  as  it 
approaches  the  shoreline.  The  amplitude  is  seen  to  grow  slowly  till  the  wave 


front  reaches  the  shoreline  and  then  the  growth  is  much  more  rapid.  Less  obvious 
in  this  figure  is  the  decrease  in  wavelength  due  to  the  increasing  local  sound 
speed.  Some  steady  state  wave  profiles  for  this  case  are  shown  in  figure  5b. 

The  arrows  attached  to  the  curves  indicate  the  direction  of  the  water  velocity 
at  the  corresponding  locations.  (The  scale  on  these  figures  is  such  that  the 
bottom  profile  has  slope  10^  and  thus  appears  almost  vertical. ) 

Figure  6a  and  6b  show  corresponding  wave  profiles  for  uj  =  7  and  a  =  10"^. 
The  shortening  of  the  wave  length  is  clearly  illustrated  here.  The  local  wave¬ 
length,  L,  at  a  position  x  is,  in  the  linear  theory, 

3.  «*>- Is  ■=«■§*  M 1/2 

As  a  rough  check  we  examine  the  curve  in  figure  6a  for  t  =  I.765  (or  cyle  No.  500). 
Calling  the  distance  between  two  successive  nodes  or  two  local  extrema  a  half 
wavelength,  L/2,  and  comparing  this  to  the  theoretical  results  (3)  evaluated 
at  the  midpoint  of  the  two  points  in  question  we  get: 


Somewhat  better  agreement  can  be  obtained  by  replacing  (3)  with  average  local 


wavelengths  defined  in 


,.a*i 


dx1  n 
c(x' )  =  u 
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a  = 

10'3  1 

max 

'  min 

a  =  2  x  10 


4.0 ! 

j 

7.0/7.9 

-7.1/-7.7 

4.5 

7.69 

-7-77 

5.0 

7-87 

-8.07 

5.5  ; 

i 

7. 0/8. 2 

-8.4 

6.0  1 

8.21 

-9.45 

6.5  i 

25/125 

j 

-22/-70 

7.0  j 

I 

11.8 

-8.5 

1 

7.5  i 

12.4 

-9.3 

8.0  j 

!  12.9/14.2 

-9.1/10.4 

I 

9.0  j 

1 

5.4/44 

-14/ -47 

.0.0 

43/80 

-18/-43 

7.70 


-7.00  7.18  -7.13  7.V7.9  5  -6.6 


-8.15  11. 2/18  ‘  -7/-IO.5I  8.2/30  -3.V-18.7 


-8.1/ -9. 3*’  5.8/33  -4.2/-221 


50/53  '  -28/-30 


Table  I 


Maximum  and  minimum  shoreline  excursions  in  "steady"  state  |i.e. 
min  for  T  "large"  and  P  at  least  a  period  . 


max 

T<t<T4P 


iLll 

a 


,  min  L  for  T  "large"  and  P  at  least  a  period 

ana  Kt<T+P  a 


6.  The  breaking  of  waves 

The  theory  of  wave  amplification  reviewed  in  section  2  concerns  the 
case  in  which  the  waves  do  not  break  and  bores  do  not  form.  But  the  excluded 
case  is  the  most  important  one  in  the  study  of  tsunamis,  since  experience  shows 
that  they  always  form  bores.  Therefore  we  shall  now  present  a  theory  of  wave 
breaking  with  the  consequent  formation  and  growth  of  a  bore.  We  shall  base  our 
analysis  on  the  nonlinear  shallow  water  theory,  assuming  the  motion  to  be  in  a 
plane  and  the  bottom  to  slope  uniformly.  The  main  mathematical  technique  to 
be  used  is  the  theory  of  characteristics,  which  will  be  combined  with  the  theory 
of  weak  shock  waves. 

A  qualitative  description  of  wave  breaking  based  on  the  nonlinear 
shallow  water  theory  was  given  by  H.  Jeffreys  in  193^«  It  is  that  the  higher 
parts  of  a  wave  travel  faster  than  the  lower  parts  so  that  the  wave  steepens  as 
it  propagates  and  ultimately  breaks.  This  explanation  was  repeated  by  J.  J.  Stoker 
[3]  and  illustrated  by  the  calculation  of  a  wave  profile  at  successive  instants, 
showing  how  it  steepens.  Stoker  also  calculated  the  time  and  place  at  which  the 
profile  becomes  vertical  at  some  point,  which  is  the  time  and  place  at  which 
breaking  may  be  said  to  occur.  These  calculations  were  for  waves  in  water  of 
constant  depth. 

In  this  section  we  shall  also  calculate  the  wave  profile  at  successive 
instants,  and  determine  the  time  and  place  at  which  breaking  occurs,  for  waves 
on  a  uniformly  sloping  beach.  In  the  next  section  we  shall  determine  the  manner 
in  which  the  height  of  the  resulting  bore  grows  as  it  propagates  toward  shore. 
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We  begin  by  vriting  the  equations  of  the  nonlinear  shallow  water  theory 
in  characteristic  form  as  follows 

1.  (u  +  2c)t  +  (u  +  c)(u  +  2c)x  -  ghx  =  0 

2.  (u  -  2c)t  +  (u  -  c)(u  -  2c)x  -  ghx  =  0 

For  a  uniformly  sloping  bottom  h  is  constant.  Therefore  we  write  gh  =  -m  where 

X  X 

m  is  a  positive  constant.  Now  (l)  and  (2)  can  be  integrated  with  the  result 


(see 

[3  ],  P.  294) 

3. 

u  +  2c  +  mt  =  k+  =  const. 

on  C+  defined  by 

SIS' 

ii 

4. 

u  -  2c  +  mt  =  k  =  const. 

on  C  defined  by 

dx 

dt  = 

The 

curves  C+  and  C  are  called 

characteristics. 

From  (3)  and  (4)  we  find 

5.  u  =  |f(k+  +  kj  -  wt 

6.  c  =  |  (k+  -  k_) 

Now  the  equations  for  the  characteristics  can  be  written 

7*  S  =  i  k+ + 1 k.  - mt  on  c+ 

8*  if =  t  k+ +  i k.  ■ mt  on  c. 

The  constants  k+  and  k_  can  be  detemined  from  (3)  and  (4)  if  the  values 
of  u  and  c  are  known  at  some  point  on  each  characteristic.  For  example  at  the 
shoreline  x  =  £(t),  we  know  that  c  =  0  and  u  =  d£/dt.  Therefore  if  a  character¬ 
istic  touches  the  shoreline  at  t  =  t  ,  we  have  on  it 

8 
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^(t) 

9.  k  =  mt  + 

+  s  dt 

If  the  shoreline  is  at  rest,  (9)  becomes 

10.  k,  =  mt 

+  s 

In  a  region  of  the  water  which  is  undisturbed  we  have  u  =  0  and 
c  =  \/-mx  .  Therefore  the  undisturbed  characteristics  satisfy  the  equation 

11.  ~  =  +  -./-rax  on  C, 

dt  —  '  + 

The  solutions of  (11)  are 

12.  t  -  t  =  +  -[(-mx  ^^-(-mx)1/2]  on  C 

o  -  m[ v  o'  J  + 

Let  us  now  consider  a  C+  characteristic  which  crosses  a  family  of  C 
characteristics  all  of  which  have  come  from  the  undisturbed  shoreline  £(t)  =  0. 

Then  from  (10),  k  =  mt  on  them  so  (5)-(7)  become 

•  s 

13.  u  =  ^(k+  +  mtg)  -  mt 

Ik.  c  =  |{k+  -  mts) 

g  -  i  K  x  t 

We  now  make  the  basic  assumption  of  the  theory  of  weak  waves,  which  is 

that  the  C_  characteristics  are  undisturbed.  Then  their  equations  are  given  by 

(12)  with  the  lower  sign.  Since  they  touch  the  undisturbed  shoreline,  we  may 

set  tQ  =  tg  and  xq  =  0  in  (12).  In  this  way  we  get  tg  in  terms  of  x  and  t.  Then 

by  using  this  expression  for  t  we  can  rewrite  (13) -(15)  as 

s 
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16.  u  =  |  k+  ‘  ‘  ('nix)1/2 

il-  c  =  \  k+  -  IT  + 

18-  8  =  i  <k+  - mt)  -  |(-mx)1/2 

To  solve  (18)  we  introduce  t  and  z(t)  defined  by 

l/p 

19.  T  =  mt  -  k^,  z(t)  =  (-mx) 

Then  (18)  beccxnes 


Since  (20)  is  a  homogeneous  equation,  it  can  be  solved  by  introducing  u(t) 
defined  by 

21.  u(t)  =  t’1  z(t) 

In  terms  of  u,  (20)  becomes 


22. 


du  _ 
dr  8u 


u 


The  variables  are  separated  in  (22)  so  it  can  be  integrated  to  yield 

.1  ,  3/5  ,  2/5 

23.  T  1  =  A1(u  -  J)  (u  +  |) 

Here  is  an  arbitrary  constant. 

Upon  eliminating  u  by  means  of  (21)  and  simplifying,  we  can  write  (23)  as 

3  T  2 

24.  (z  -  f)  (z  +  §)  =  A 


Here  A  i.6  another  constant.  In  terms  of  t  and  x,  (24)  becomes 


IK 


(-mx) 


1/2  + 


St  k.f 

2  '  2  J  _ 


25. 
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The  constant  A  can  be  expressed  in  terms  of  some  point  x  .  t  on  the  char  a  c  ter - 

o  o 

istic,  and  ^25)  can  be  rewritten  as 


The  constant  k+  is  given  by  (3)  in  terms  of  the  values  of  u  and  c  at  (xq,  tQ)  as 

27.  k+  =  u(xQ,to)  +  2c(xo,to)  +  mtQ 

Thus  u  and  c  at  any  point  (x,t)  satisfying  (26)  are  given  by  (16)  and  (17), 

with  k.  given  by  (27). 

+ 

To  solve  (26)  for  x  in  terms  of  t,  we  first  define  )  and 

f(xo’to)  by 

28.  c(xo,to)  =  (-tnxo)1/2  +  f(xo,tQ) 

29.  f(x0>t0)  =  u(x0>to)  +  2fl(xo»to) 

Now  if  f(xQ,tQ)  =  0  the  squared  term  on  the  right  side  of  (26)  vanishes.  By 
equating  to  zero  the  squared  term  on  the  left  side  of  (26)  we  rederive  the 
equation  (12)  for  an  undisturbed  C+  characteristic. 

Now  let  us  suppose  that  f(xQ,tQ)  is  not  necessarily  zero,  but  is  small. 
Then  we  may  expect  the  characteristic  (26)  to  differ  only  slightly  from  the 
undisturbed  characteristic  (12).  Therefore  we  shall  try  to  express  it  in  the 
form 


30. 


-  to  =  |[(-mxo)l/2  "  ("mx)1^2J  +  4 
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Here  q  is  to  be  determined  by  substituting  (30)  into  (26).  By  keeping  terms  of 
the  zeroth  and  first  powers  in  f(xo,tQ)  we  find  that 


Thus  for  small  motions,  the  C+  characteristic  is  given  by 
32.  t-t 


m« 


r|X>3A  i 

]-;Vo> 

L\*  J 

We  may  now  use  (32)  in  (l6)  and  (17)  to  obtain  u  and  c  on  the  C+ 
characteristic  (32),  with  the  results 

X  3/ 

33.  u(x,t)  =  |  f(xo,t0)  ^ 

X  3/4 

3lt.  c(x,t)  =  (-mx)1'2  +  |  f(xo,to)  ^ 

These  results  are  not  valid  near  x  =  0  because  there  it  is  not  correct  to  omit 
higher  order  terms,  as  was  done  in  deriving  (31)*  In  fact  the  full  expression 
(26)  for  the  C+  characteristic  leads  to  finite  values  of  u  and  c  everywhere 
including  x  =  0. 

It  is  easy  to  see  from  (32) -(3^)  that  wave  crests  steepen  in  front  anu 
become  vertical  as  they  approach  shore.  To  show  this,  let  us  evaluate  u  (x,t) 
assuming  that  the  initial  values  of  u  and  c  are  given  on  the  line  xq  =  constant. 
Then  from  (33)  we  have 


35.  ux(x,t)  =  -  g  f(xo,to)x3/4  x’l/k  +  \  ft  (x0,t0)(x0/x)3/U  -g 


Now  from  (32)  we  find 
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.  2/4  -» 

36.  .  £  ,  (^,)-l/2  +  jL  f(s  >t  )x3A  x-7A  .  i  f  (,  ,t  )  -ip: 

dx  x  '  4m  v  o'  o'  o  m  t  v  o*  o  \x  dx 


LJ 


Solving  (36)  for  dtQ/dx  yields 


37-  s2  -  ;  ft  I'V^374-1]-1  1  <-“)'1/2  * 


3 iJ'k  ^ 

km 


We  see  from  (37)  that  dt^/dx  becomes  infinite  when  the  denominator  vanishes  which 
occurs  if  u  +2f.  >  0,  Then  from  (35 ),  u  also  becomes  ininite.  The  value  of 

U  L  X 

c  o 

x  at  which  this  occurs  is  found  from  (37)  to  be 


38. 


x  =  x 


l  0-> 

The  time  at  which  this  occurs  is  obtained  from  (32)  by  inserting  the  value  (38) 
for  x.  The  first  occurrence  of  an  infinite  value  of  u  is  determined  by  the  value 

X 

of  tQ  which  minimizes  t.  Exactly  the  same  type  of  consideration  applies  if  the 

initial  values  of  u  and  c  are  given  at  t  =  constant  or  on  any  other  curve  in 

the  x  ,  t  plane. 
o’  o 

When  u  becomes  infinite,  c  also  becomes  infinite.  The  first  occurrence 
x  x 

of  this  vertical  slope  for  a  given  wave  is  called  breaking  of  the  wave.  There¬ 
after  the  solution  (32) -(3^)  is  multiple  valued.  This  multiple  valuedness  can 
be  eliminated  by  the  introduction  of  a  discontinuity,  called  a  bore,  into  the 
subsequent  wave  profile.  We  shall  now  show  how  this  is  to  be  done. 
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7.  Hie  formation  and  growth  of  bores. 

A  bore  is  a  discontinity  in  '  vfcer  depth  and  velocity  across  which  two 
jump  conditions  must  be  satisfied.  These  are  the  mass  and  momentum  equations 
which  may  be  written  in  the  form 
1.  +  h)(ux  -  6)  =  (t)2  +  h)(Ug  -  b) 


2. 


(il1-rti)u1(u1-S)+  |(r)14h)2  =(ri24h)u2(u2-S)+  |(tj2-Hi)2 


Here  and  u^  denote  the  wave  height  and  water  velocity  on  side  one  of  the  bore, 
s  is  the  bore  velocity  and  u^  denote  quantities  on  side  two.  From  (l)  and 
(2)  we  find  for  a  weak  bore,  i.e.,  one  for  which  «  h. 


3- 


S  =  u  +  </ih 


3T)l_Tl2 

l  +  -fe r + 


b. 


U1  U2 


+  M 


v»2 


If  the  fluid  crosses  the  bore  from  side  two  to  side  one,  the  plus  sign  applies 
in  (3)  and  (4),  and  side  two  is  called  the  front. 

When  side  two  is  the  front,  we  find  from  the  definition  of  c  and  from 

(4)  that 


5. 


Cl+Ul+C2+U2  =  2 


+ 


Vsfc  1 


i  + 


3\"T>2 


Upon  comparing  (5)  with  (3)  in  which  the  plus  sign  is  chosen,  we  find  that 
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6. 


S  =  |  (c^-K^+Cg+Ug)  +  .. 


2  \  dt  dt  j  ’ 


Thus  the  speed  of  a  weak  bore  is  the  average  of  the  C+  characteristic  speeds 
on  its  two  sides.  This  result  also  holds  in  gas  dynamics  and  in  certain  other 
cases. 

Hie  expression  (6)  for  the  velocity  of  the  bore,  together  with  the 
expression  (6.32)  for  the  C+  characteristics  suffices  for  the  determination  of 
the  location  of  the  bore  as  a  function  of  time.  To  determine  the  bore  location 
we  again  use  the  notation  z  =  (-mx)1^  and  we  also  introduce  c  =  .  We 

assume  that  the  incident  wave  is  given  at  x  =  x  which  corresponds  to  z  =  z  . 

o  o 

Suppose  that  two  C+  characteristics  which  leave  x  =  xq  at  t  =  t^  and  t  =  tv> 
with  t^  >  t^,  meet  at  the  bore.  The  equations  of  these  characteristics  are 
given  by  (6.32)  with  t  replaced  by  t  or  t  .  In  terms  of  z  and  0  these  equations 

are 

7.  a  +  z-zQ  =  a1-|  f(o1)  [(zo/z)3/2  -  l] 

8.  0  +  z  -  zo  =  o2  -  |  f(a2)|(zo/z)3/2  -  l] 

Here  f(a^  denotes  f(xQ,t1). 

Let  us  write  the  equation  of  the  bore  in  the  form 

9.  a  +  z  -  zq  =  G(z) 

To  determine  G(z)  we  first  combine  (7)  and  (8)  with  (9)  to  obtain 

10.  G(z)  =  Oj.  _  2  fK)[(*0A)3/2-l! 

G(z)  =  °2  ~  2  f(a2'  O^c/2^2"1] 


11. 


-25- 


Next  ve  wish  to  utilize  (6).  Before  doing  so  we  observe,  from  the  definitions  of 
z  and  o ,  that 


.  dx  2z  dz  dz 

S  ”  dt  ”  "  m  dt  “  ~z  da 


Now  from  (9) 


dz _ 1 

dO  n  \  t 


uu  G '  ( z )  -1 

Thus  (12)  and  (13)  yield 


-1  -  G'(z) 


s  =  z[l  +  G*  (z)] 


To  obtain  dx^/dt  and  dx^/dt  we  differentiate  (7)  and  (8)  to  find 

»•  £  -  ^  ■  -W*  i  .*V.  .u*  i  «,,>*, f  S'h 

!6.  i  S'2)-Kziu  \  t  S'2) 

Upon  using  (l4)-(l6)  in  (6)  we  obtain 

17.  2G'  (z)  =  |  zZJ2  z'5/2ff(a1)+f(o2)] 

The  problem  is  now  that  of  solving  (10), (11)  and  (17)  for  o^,  op  and 
G  as  functions  of  z.  To  solve  this  problem  we  follow  the  method  of  G.  B.  Witham 
[4],  We  differentiate  (10)  and  (ll)  with  respect  to  z  and  add  the  results, 
obtaining 

18.  2G'  (z)  =  J  Z3/2  z‘5/2ff(a1)+f(a2)J  +  ^  {1-  |  f '  (a^  [(zo/z)3/2  -  1 j  ) 


+  d f  (l-|f’(a2)[(Zo/z)^"  -1]} 
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By  taking  account  of  (17),  we  can  simplify  (18)  to 


3T  ■, 1  -  -T2- 


[<«>)3/2-^  ?  -1-  tv.)3/2-]  ■ 


To  eliminate  z  from  (19)  we  subtract  (11)  from  (10)  and  find 


(«  A)3/2  - 1 . 


Now  (19)  can  be  rewritten  as 


0,  -  a 


V*u  1  -  Urt  uu  _  (J.  “  O- 

dT  ^-f'(°i)  f(^)-f(c2)  +5T  (  1  '  f,(o2)  iT^ff 


It  is  convenient  to  multiply  (21)  by  f(o^)-f(a2)  and  then  to  earrange 


it  in  the  form 


d°o  daT  T  A 

f(c2)  d T  -  f(ol}  dT  =  2  h  {[°2-o1][f(a2)+f(a1)]) 


Integration  of  (22)  yields 


c. 

f  f(a)do  =  |  (^2”a1)  [f(a2)+f(a1)] 


The  constant  of  integration  in  (23)  has  been  chosen  so  tha+  ^23)  is  satisfied 
by  0^  =  o^,  which  corresponds  to  the  first  point  on  the  path  of  the  bore.  The 
condition  (23)  determines  the  correspondence  between  and  0^*  The  geometrical 
interpretation  of  (23)  is  that  the  area  bounded  by  the  curve  f(a),  the  a-axis  and 
the  lines  a  ==  and  a  =  <?2  is  equal  to  the  area  of  the  trapezoid  with  vertices 
at  [oy  f(a1)],  [p2,  f(a2)],  [o^  0]  and  [a2,  0] .  Alternatively  we  may  say 
that  the  signed  area  between  the  curve  f(a)  and  the  chord  Ircm  [0^,  ffa^)]  to 
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[a2,  f(a2)]  is  zero*  Thus  the  chord  must  cut  off  equal  areas  on  the  two  sides 
of  the  curve. 

Once  cr2  is  determined  as  a  function  of  from  (23),  the  path  of  the 
bore  is  found  by  eliminating  G(z)  from  (9)  and  (10)  to  obtain 


24.  0-0-,+  z-  zQ  =  -  |  f(c1)  [(zQ/z)3/2  -  l] 

Now  z  is  given  in  terms  of  and  by  (20).  Thus  (20)  and  (24)  are  parametric 

equations  for  the  path  of  the  bore  with  a ^  as  parameter.  These  equations  may 
be  written  more  explicitly  in  the  form 


25. 


z  =  z 


2  (a1  - 
1  +  f(a1)-f( 


a2} 


°2) 
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°1  ‘  a2 

26.  a  =  a1+  zo  -  z  -  f(o1) 

The  discontinuities  in  u  and  e  can  be  found  from  (6.33)  and  (6.32). 

They  are 

Z  3/2 

27'  °ru2  ■  I  r, 

,  z  3/2 

28.  Cj-C2  =  2;  T-  [f(o1)-f(o2)] 

By  using  (20)  we  may  rewrite  these  results  as 

29.  ux-  u2  =  cx-  c2  =  0^-  02+  |[f(a1)-f(o2)] 

The  equation  (23)  can  be  solved  readily  for  in  terms  of  if  f(a) 


is  an  odd  function,  i.e.,  if 
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30.  f(-a)  =  -f(a) 

It  is  then  easy  to  see  that  the  solution  of  (23)  is 

31.  a2  =  -o1 


Then  (25) >  (26)  and  (29)  become 

2a,  1‘2/3 


32. 


z  =  z 


1  4 


33.  a  =  zQ  -  z 

j4.  ul'u2  =  ci‘c2  =  2ol  +  f^°l^ 

To  illustrate  the  results  (32) -(34)  let  us  apply  them  to  a  sinusoidal 
incident  wave  for  which 

35.  f(a)  =  A  sin  ^ 


Then  (32)  and  (34)  become 


36. 


z  =  z 


2o, 


1-2/3 


1  + 


A  sin 


2nol 


2x1^ 

37.  ux-  u2  =  cx-  c2  =  20l  +  A  sin  - 

We  can  simplify  these  formulas  by  introducing  in  place  of  the  parameter  6 
defined  by 

2na, 

38.  &  =  t 


Then  (36)  and  (37)  become 
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39. 


z  =  z 


_s_r2/3 

sin  5  J 


40. 


Hr 

u,  -u_  =  C-.  -c_  =  A  —  5  +  sin  5 
12  12  |_xtA 


T  r2li 


From  (39)  ve  find  that  the  bore  forms  at  z  =  z  (1+  =rr)  where  5  =  0 

O  xtA 

and  travels  toward  the  shore  z  =  0  which  it  reaches  when  5  =  xt.  From  (33)  we  see 

T 

that  the  bore  path  is  an  undisturbed  C,  characteristic  starting  at  o  =  z  -z  (1+  — ) 

+  c  ov  xtA 

and  reaching  the  shore  at  a  =  zq.  The  maximum  discontinuity  across  the  bore 

occurs  at  5  =  cos  ^(-T/xta)  provided  T/xtA  <  1,  and  its  value  is 

ol 

I  m  “i  m  I  ;  m 

41.  ui'u2  =  cl’c2 


f  T  -1 

.  -T 

/  w  2 

1_ 

\^C0S 

IrtA/ 

l 

If  T/xtA  >  1,  the  maximum  discontinuity  occurs  at  6  =  xr  i.e.,  at  the  undisturbed 

gT  tan 

shoreline,  and  its  magnitude  is  T.  In  terms  of  the  original  variables  T  =  - - — 

where  a  is  the  beach  angle  and  Tq  is  the  period  of  the  incident  wave. 

When  T  >  xtA  the  ratio  of  the  maximum  bore  discontinuity  to  the  incident 
amplitude  is  T/A  =  gTQ  tan  a/2A.  Thus  this  factor,  which  we  may  call  the  velocity 
amplification,  increases  with  the  period  of  the  wave  and  with  the  beach  angle. 
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